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DYNAMICS OF TWO SLOWLY ROTATING POINT -MASS VEHICLES 
CONNECTED BY A MASSLESS TETHER AND 
IN A CIRCULAR ORBIT 

By William M, Adams, Jr, 

Langley Research Center 

SUMMARY 

A study has been made of the dynamics of two point masses in circular orbit con- 
nected by a massless tether and rotating relative to the local vertical with angular rates 
of the same order of magnitude as the orbital angular velocity. Only motion in the orbital 
plane is considered. Equations of motion are developed which specify the behavior of the 
system during periods when the tether is in tension and slack. Results that allow predic- 
tion of the occurrence of a slack tether are obtained. It is also shown that the addition of 
damping, coupled with oscillations between a slack and taut tether, causes the system to 
seek closer alinement with the local vertical, 

INTRODUCTION 

Use has been made of a tether in Gemini flights for extravehicular maneuvers and 
as a means of providing artificial gravity. The tethered configuration (fig, 1) has also 
been proposed as a means of station keeping. (See ref, I*) The characteristics of the 
motion of tethered bodies have been the subject of several theoretical studies in recent 
years, (See refs, 2 and 3.) These studies considered systems that were spinning with 
angular rates that were much larger than the orbital angular velocity. 

If the system is spinning relative to the local vertical with an angular rate that is of 
the same order of magnitude as the orbital angular velocity, the gravity gradient between 
the two masses has a pronounced effect upon the dynamics of the system. This study is 
an analysis of the behavior of such a configuration* 

Several simplifying assumptions are made in carrying out the study. The vehicles 
are taken to be point masses constrained to the orbital plane of the center of mass which 
is assumed to be in a circular orbit about the earth. Additional assumptions are that the 
tether is massless and can be characterized by a force constant, and any damping present 
in the system is viscous. 



Equations of motions are developed which specify the behavior of the system during 
periods when the tether is in tension and slack. Results which allow prediction of the 
occurrence of slackness in the tether are presented. The study also shows that if damping 
is present and oscillations between a slack and taut tether occur, the motion of the masses 
relative to their center of mass is affected in such a way that the system becomes more 
nearly alined with the local vertical and no further slackness occurs . This result is 
desirable if gravity- gradient stabilization is an objective. 

SYMBOLS 


D viscous damping coefficient 

F elliptic integral of first kind 

h unit vector normal to orbital plane in direction of orbital angular momentum 

I moment of inertia of taut system about its center of mass, + n^^ 

i,j,k right-handed set of unit vectors in direction of local horizontal, local vertical, 

and normal to orbital plane of center of mass, respectively (see fig. 4} 

K constant depending upon initial conditions (defined in eq. (8)) 

k spring constant of tether 

L Lagrangian 

l unstretched length of tether 

unstretched lengths of tether from center of mass to mj and mg 

mj,m 2 masses at end of tether 

n = r - l 

P period of oscillation of a simple spring assuming viscous damping 

R vector from center of earth to center of mass 
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> 1*2 

rj,r 2 

F= in- 


vertors from center of earth to 
vectors from center of mass to 

F 2 



and 

m 2 

m l 

and 

m 2 


T 

t 


K>v<p 



V 1> V 2 

F 1> F 2 


tension in tether 
time 

unit vectors in radial and tangential directions, respectively 

inertial velocity of center of mass 

inertial velocities of m^ and m 2 

velocities of mj and m 2 relative to center of mass 


v = V-i - 


1 ' v 2 


x,y,z 


components in i,j,k directions 


fraction of radial velocity remaining after a slack-taut- slack transition 


$ 

<P 


angle between ref erence line in orbital plane and local vertical 


product of universal gravitational constant and mass of earth 


angle between local vertical and r 




angle defined as sin" Y - sin (p\ 

fK+| 

2 


Subscripts: 

B bound upon angular oscillation 

£ at an extremum 
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I 


with respect to inertial space 


s value of quantity when tether becomes slack 

T value of a quantity when tether reaches full extension after having been slack 

o at initial time 

1,2 refers to masses and respectively 

r ,<p in radial and tangential directions, respectively 

Dots over symbols denote derivatives with respect to time, A bar over a symbol 
denotes a vector whereas a caret over a symbol denotes a unit vector. When a vector 
quantity is written without a bar, the magnitude of the quantity is denoted, 

THEORY AND ANALYSIS 

Sets of equations are developed that specify the motion of the tethered system during 
periods when the tether is slack and taut, and a method for predicting whether an initially 
taut tether will become slack is derived. 

Dynamics of the System During Taut Tether Periods 

Consider the system shown in figure 1, Two point masses are connected by a mass- 
less tether and their center of mass is moving in a circular orbit with constant orbital 
angular velocity of magnitude 9. The tethered system is rotating (in the orbital plane) 
relative to the local vertical with angular rate 

The equations of motion of the system are obtained by using the Lagrangian approach 
with the separation of the two masses r and the angle from the local vertical to the line 
connecting the masses cp serving as generalized coordinates. The Lagrangian is defined 
as 

L = Kinetic energy - Potential energy 
and the equations of motion are given by 



and 
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d_/&L\ 3L_ Di 

dtUJ' 9r 


where D is a viscous damping coefficient that is assumed to characterize the damping 
present in the system. 


The Lagrangian may be written as 


1 ■ l m Wi • W, + 1 “N • W, * ^ * it - 1 k(r - ,)! (rS0 (3) 

In this expression and {^ 2 )^ are the velocities of nij and m 2 relative to 

inertial space, Rj and R2 are the distances to mj and m2 from the center of the 
earth, l is the unstretched length of the tether, M is the universal gravitational con- 
stant times the mass of the earth, and k is the tether force constant. The quantities 
Rf, R 2 , V lf and V2 may be written in terms of R, r, V, v, l, and <p as 


R 1 = Jr 2 + rj 2 + 2r jR cos cp = JR 2 + r 2 + 2^rR cos ® 


R 2 = ^R 2 + r 2 2 - 2r 2 R cos <p ~ /r 2 + j r 2 - 2f-^jrR cos <p 


Vj « Oh x ^R + r + Vj = ^(?R sin <p + rjtTr + cos <p + -j- r(e + <p) 1 


V 2 = j9h x + r 2 ^ + v 2 - ^£7R sin 0 - -j- r^U r + 6 


dR cos cp - —■ r 


• 1 . 

(fl + 0 ) 1 


In these equations, 0h is the constant orbital angular velocity of the center of mass, R 
is the radius from the center of the earth to the center of mass, U r is a unit vector from 
m 2 toward and II ^ is a unit vector defined as - h x U r . Substitution of equa- 

tion (3) into equations (1) and (2) yields the equations of motion of the tethered system: 


r%(p 4- 2rr(f? + ip) + ^ r 2 # 2 sin 20 = 0 


2 * & 
r + D ~ r + r0 2 - 3fS 2 r cos 2 <p - r(fi + <p) 2 + ~-(r - 1) = 0 (rl?) 


•\2 , kl 2 
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where I = mjfj 2 + The assumptions made in deriving these equations were 

the following; 

(1) The center of mass remains in a stationary circular orbit 

(2) The massless tether can be represented by a force constant, and any damping 
present is viscous 

(3) In the expansion of ^R^Rj and , terms involving (r/R) raised to 

powers higher than 1 can be neglected. 

The tension in the tether produced by the gravity gradient and the slow angular rates 
under study corresponds to extension of the tether that is insignificant relative to the 
tether length, and the tether can be regarded as of constant length. Under this condition, 
equations (4) and (5) simplify further to 

ip~ ~ ~ 0 2 sin 2 (p (6a) 

a 

T = k(r - l) = i|(fl + <f>) 2 - 0 2 (l - 3 cos2<p| (TiO) (6b) 

This simplification is not valid, in general, immediately following the period when the 
tether again becomes fully extended after having been slack; in such cases the relative 
velocity between the two masses can produce tension that is considerably higher than that 
produced by the gravity gradient and the slow rotational rate. Analysis of the stretching 
motion during such a period will be presented in a later section of the paper. 

Equation (6a) specifies the angular acceleration of the tethered system about its 
center of mass produced by torque resulting from the gravity gradient between the two 
masses. Solution of this equation specifies the angular orientation relative to the local 
vertical at arbitrary times. The expression 

ip = - ~ 0 2 sin 2(p 


may be integrated once to yield 



(7) 


where 

clearly 


K is a constant depending upon the initial angular velocity and orientation and 
obeys the inequality 


K = 



~ cos 2 5 - 
2 o 


3 

2 


( 8 ) 
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By disregarding for the moment the possibility that the tether may lose tension 
and examining equations (6a) and (7), several statements can be made regarding the rigid- 
body motion. These statements are summarized as a function of the parameter K in 
table I. 


TABLE I.- TYPES OF MOTION WHEN TENSION MAINTAINED 

Type of motion 

System remains alined with local vertical 
System oscillates relative to local vertical with amplitudes 

System approaches local horizontal and requires Infinite time 
to get there 

System makes complete revolutions relative to local vertical 


It should be emphasized that table I is applicable for the tethered system only while ten- 
sion is maintained and must be modified somewhat to account for the occurrence of a 
slack tether. 

Equation (7) may be integrated in closed form in terms of elliptic integrals of the 
first kind. The exact form of the result will depend upon the value of K- 

The closed-form expressions are derived for the following cases: 

Case I.- K = - 1 {<p = 0,jt; <p = ip = 0) 

£t 

Case H.- -|< K<| 

6 b 

Case in.- K = | 

2 

Case IV. - K > | 

3 

For case I for which K = - inspection of equation (6a) reveals that the system 

& 

tends to maintain this orientation because any small angular rate produced by an external 
source will be opposed by the torque created by the gravity gradient with resulting oscil- 
lations about the local vertical as seen in case EL 



Value of K 
3 

"2 


2 <K< 2 
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For case E where - | < K < |> by substituting the relation 



into equation (7) one obtains for the time elapsed during the rotation of the system from 





K + 


s «w - 




where F is an elliptic integral of the first kind. This equation reveals that oscillations 
relative to the local vertical occur provided the tether is taut. The amplitude of the oscil- 
lations is 




<sin-l(i /^" 2 


For ease El where K = |, the integral becomes 

At = ^ = p[ F(1,w ' F ( 1,,? °3 


#? I08e (i 


tan (p + sec <p \ 
tan <p Q + sec <pl j 


In this case an infinite amount of time would be required for the system to rotate from a 
particular orientation to <p = r/2 or -ir/2. Also it can be seen from equation (6a) that if 
the system were at rest at <p = ±tt/2, it would tend to diverge from this point since any 
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angular rate relative to the local vertical created by an external torque would lead to com- 
plete revolutions of the bodies about the local vertical as is shown in the following case. 

For case IV, K > As can be seen from equation (7), rigid body motion is one of 
complete revolutions relative to the local vertical. The integral of equation (7) is 



Since each of these integrals expresses the orientation of the tethered system as a 
function of time only if the tether is taut, it is necessary to examine the expression for 
tension to determine when the tether may become slack. Substitution of equation (7) into 
equation (6b) yields the following expression for the tension in the tether: 


T = j 0^(k + ^ - 3 sin2<p ± 2 + j| - 3 sir&cp + 3 cos2<pJ 

By finding the angular orientation at which T is minimum as a function of K and 
evaluating T at that point, one can determine whether slackness will occur. The 
extrema can be found by solving for the roots of — . Forming this derivative gives 

dip 


(9) 


dT I *2 

~~ = -6 ~ 9 sin cp cos <p 

dip l 


2 ± + ~ - 3 sin2<pj 


- 1/2 


(10) 


Therefore ^- = 0 if <p = 0, tt/ 2, n, 3ir/2. In addition = 0 if (p < 0 provided 


dtp 


dip 


-( K + !- 3 si " Ss ) 


- 1/2 


= 0 


The latter condition can only be satisfied for a certain range of initial conditions, that is, 
for a certain range of K. This statement is true because of the restriction that 

| sin ip ^ I = ^ s 1 and real. These restrictions can be written in terms of K as 

-!= k 4 

4 4 
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where the positive limit requires Jsm <p^,j ^ 1 and the negative limit requires that 

sin be real, 

E 

It is necessary to determine whether these extrema are maxima, minima, or 

2 2 
neither. If evaluated at <p = ©„ is negative, T is maximum; if i S posi- 

dcp* dtp 2 

tive at (p = <p^, T is minimum. Upon applying this test, the results given in table II 
are obtained. 


TABLE II. - ORIENTATIONS AT WHICH TENSION IS EXTREMAL 





Type of extremum 

sin"l 



T 

Is minimized 

0, jr; <p E > 0 

T 

is maximized 

0, tt; cp E < 0; - | < K 

T 

is maximized 

0, tt; <p E < 0; - | < K < - | 

T 

is minimized 



„ > 0; K > - 

T 

is minimized 


2’ 2 

E ’ 2 




£ 3£- a 
2’ 2 ’ ' 


T 

is minimized 


£ A, 

2’ 2 ’ 

< 0; | < K < 1 

T 

is maximized 


Figures 2(a) and 2(b) present the normalized tension as a function of the angle from 
the local vertical for several values of the parameter K. It is seen that the tether 
remains in tension if K<l,if K > ^ for <p positive, or if K > U f or negative. 

Note the consistency between figures 2(a) and 2(b) and table n as to the type of extrema. 
When tension is negative, the curves are dotted to emphasize that such portions of the 
curves do not represent the dynamics of the tethered system. 

The angle at which the tension becomes zero as a function of K is shown in fig- 
ures 3(a) and 3(b) for <p Q < 0 and > 0, respectively. The results from this figure 
can be used to modify table I as shown in table HI. 
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TABLE EEI.- MODIFICATION OF TABLE I TO DENOTE CONDITIONS 
FOR OCCURRENCE OF SLACKNESS 


Value of K 

Type of motion (limited to orbital plane) 

Ks= “i 

System remains alined with local vertical and in tension 


System oscillates relative to local vertical with tension 

-§<KS1 

( K+ l\ //K + |\ 

being maintained, -sin" 1 ll / — J 1 cp S sin' My — ^ — J 

1 < K “ | 

Slackness occurs in tether 

K > |, (f> > 0 

Complete revolutions relative to local vertical occur 
with tension being maintained 

|<K<^s i><0 

Slackness occurs in tether 

K>1±,<P<0 

Complete revolutions occur relative to local vertical 

with tension being maintained 


Slack Tether Period 


If the tether is initially slack or becomes slack during rotation, the two masses will 
(again neglecting the effect of the tether) move independently of each other as long as their 
separation is less than the length of the tether. An outline of the development of the equa- 
tions of motion as given in appendixes A and B of reference 4 is presented. 

Consider the system as shown in figure 4. The Lagrangian associated with is 

given by 

<u) 

The x- component of the equations of motion is 



with identical equations for y^ and Zj. Performing the indicated operations and making 
use of the property that the orbital angular momentum is constant yields the equations of 
motion relative to the center of mass: 
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•2 „■ • •• M x i 

xj - Xjfl - 2y x 9 - y x B + — | = 0 


(X >Kh + R ) 


y 1 -y 1 »‘ + 2x 1 9+» 1 8-| !r +- l . ‘ ■ 0 


i 1+ £4- 0 


Specializing to the case for which the two masses are in the plane of the circular orbit of 

the center of mass, expanding Rj = jxj 2 + in a binomial series, and 

retaining only first-order terms in x^ and yj results in the following set of differen- 
tial equations: 

xj - 20yj = 0 


y t + 20xj - 30 2 yi = 0 

These approximate equations may be immediately integrated to yield the position and 
velocity of mj relative to the center of mass: 



The equations specifying the position of m 2 are identical in form to these equations. 
Thus, equations depicting the motion of the two masses during the slack tether period 
have been presented. 


Behavior During Slack- Taut- Slack Transitions 

1 

After a period of time the slack tether will again become fully extended, and relative 
motion between the two bodies in the radial direction will not, in general, be zero. The 
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resulting tension in the tether can be considerably higher than that produced by the gravity 
gradient between the two masses, and the assumption that the length of the tether remains 
constant should no longer be made. One must consider, therefore, the more complete 
radial equation of motion (eq. (5)): 

r + D j- r + rg2 _ 3g2 r co s%<p - r(g + + k ^-{r - Z) = 0 (r ^ l) 


By making the substitution n = r - l, this expression becomes 

n + D ^ n + jk + (pf 1 - 0 2 (l - 3 eos2<p) |n = + cp j 2 - e?2(i - 3 cos^Jji: 

For the low angular rates considered in this study » (0 + - 9 ^( 1 - 3 cos^g?). 

Thus, 

h + n + n = ^ + <£>) 2 - 9 ^{\ - 3 cos^jj; = ^ T 


(15) 


The term to the right of the equality is the effect of the gravity gradient and the angular 
rotation upon the extensional motion of the tether. Since the effect of this term upon the 
tether length was neglected in the earlier analysis, it will also be neglected here. Then 
the equation of motion becomes a linear homogeneous differential equation with constant 
coefficients : 

.. /2 , 1 , ;2 

n + D y- n + n = 0 


and has the solution 



where n T is the relative velocity between the two masses at the time the tether becomes 
fully extended after having been slack bp. 

The amount of damping required to remove a specified fraction of the relative veloc- 
ity between the two masses during each slack-taut-slack transition may be determined 
from equation (16). Let the relations 



(0 S« S 1) 
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and 


"(f)- “T * 0 

specify the desired conditions at the end of each damping cycle where 


P B 



Solving for the ratio of damping required to critical damping yields 



(17) 


The percent of critical damping required as a function of e as determined from equa- 
tion (17) is shown in figure 5. 

When initial conditions are such that the tether becomes slack, the motion is made 
up of periods of slackness alternating with short periods of tautness when the tether again 
becomes fully extended. A typical time history of the motion during the first slack period 
is presented in figure 6 for a case in which the center of mass is in synchronous orbit 
about the earth. It is evident from figure 6 that the relative velocity between the two 
masses is not zero when the tether again becomes fully extended and this condition results 
in a short period of tautness followed by another period of slackness. 

After a series of alternations, with damping in the system, the relative velocity 
between the two masses at full extension essentially vanishes and the system again behaves 
as a rigid body. The number of alternations before this condition is reached can vary 
from 5 or 10 to more than a hundred depending upon the damping present and the initial 
conditions. The subsequent rigid- body motion can then persist indefinitely or slackness 
can recur depending upon the angular rate and orientation of the system when it becomes 
taut. 


A computer program was written that describes the dynamics of the system in order 
to determine whether the final motion would be bounded rigid-body oscillations relative to 
the local vertical. Two tethered 930-kilogram masses were taken to be in a circular syn- 
chronous orbit, and solutions were obtained that covered the entire range of initial condi- 
tions for which slackness would occur. Time histories were computed with K 0 , the 
parameter which defines the initial conditions, chosen at increments of 0.025. Other fac- 
tors varied were damping (2 to 7 percent of critical damping), tether length (2500 and 
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10 000 meters), and the extensional rigidity (333 S kf = 5540 newtons) of the tether. 

Cases in which slackness did not occur and tumbling continued were not studied. 

The results of this analysis are shown in figures 7(a), 7(b), and 7(c). The ordinate 
is the final amplitude of oscillation of the system relative to the local vertical; the 
abscissa K 0 defines the initial conditions when coupled with the requirement that 
< 0. By reference to figures 2(a) and 2(b), it is determined that the values of Kq 
in figure 7 correspond to large- amplitude oscillations or complete revolutions relative 
to the local vertical. The dashed line in the figures is the angular bound below which 
tension is maintained and above which slackness occurs for bounded oscillations and cor- 
responds to the curve K = 1 in figure 2(a). With the exception of the absence of curves 
near Ko = 5.5 the figures show that the final steady- state condition is one involving 
bounded oscillations and no further slackness. (Final tension is instantaneously zero for 
cases with Kq near 1.) 

Increased damping tends to smooth out the curves and generally leads to a smaller 
final amplitude of oscillation relative to the local vertical. Comparison of figures 7(a) 
and 7(b) indicates that the final amplitude of oscillation is insensitive to changes in the 
tether length over the range of variations studied. Comparison of figures 7(b) and 7(c) 
shows that the same statement is true for variation in the extensional rigidity of the 
tether. 

No points are plotted in the vicinity of Kq = 5.5 because the system remained in 
such a condition that slackness would recur indefinitely. In all other cases studied, how- 
ever, the final condition is one in which the system is closer to being alined with the local 
vertical than initially and no further slackness occurs. A more complete analysis con- 
sidering the extensional motion produced (when the tether is in tension) by the slow vari- 
ations in the gravity gradient and angular rate could show enhancement of this tendency 
toward gravity gradient stabilization. 

It is of interest to observe the separation between the point masses during the slack 
periods to determine whether collisions between finite -sized vehicles might occur. The 
separation between masses for e = 0.90 and a 10 000- meter tether having a spring con- 
stant of 0.534 N/m was observed from the point of initial slackness until the rigid body 
mode was again reached, no further slackness occurring, for a range of initial values of 
K. The results, shown in table IV, are the number of times the separation was observed 
to be less than each indicated distance for a sample rate normalized to one observation 
every 8 minutes. The separation became less than 10 percent of the tether length for a 
number of cases. However, in no instance did the separation become as small as 1 per- 
cent of the tether length. Table IV was generated under the assumption that the relative 
velocity between masses was zero at the point of initial slackness. 
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TABLE IV.- NUMBER OF TIMES SEPARATION WAS 
LESS THAN THE INDICATED DISTANCES 
H = 10 000 meters] 


K 0 

Number of times separation 
was less than v/l of - 


1.0 

0.5 

0.1 

0.01 

1.40 

203 

0 

0 

0 

1.47 

I 225 

0 

0 

0 

1.60 

274 

0 

0 

0 

; 1.70 

339 

o 

0 

0 

1.80 

: 360 

14 

0 

0 

1.90 

438 

21 

0 

0 

2.00 

! 407 

27 

0 

0 

2.10 

496 

61 

7 

0 

2.20 

541 

88 

0 

0 

2.30 

573 

97 

8 

0 

2.40 

546 

43 

0 

0 

2.50 

871 

122 

0 

0 

2.75 

535 

78 

0 

0 

3.00 

566 

82 

4 

0 

3.25 

670 

116 

0 

0 

3.50 

785 

124 

5 

0 

3.75 

729 

89 

9 

0 

4,00 

632 

127 

4 

0 

4.50 

756 

107 

8 

0 


CONCLUDING REMARKS 

Equations have been developed for a tethered system for the idealized case of a 
massless tether that allow prediction of the occurrence of slackness; furthermore, it is 
shown that the addition of damping, coupled with the occurrence of slackness, affects the 
motion of the masses relative to their center of mass in such a way that the system 
becomes more nearly alined with the local vertical and no further slackness occurs. 

In addition, the minimum separation between masses was greater than 1 percent of 
the tether length for all cases having zero relative velocity between masses at the 
point of initial slackness. A more complete analysis considering the slight exten- 
sional motion produced by the variation in the gradient and angular rate occurring 
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while the system is in tension could show enhancement of the tendency toward gravity- 
gradient stabilization. 

It should be pointed out that the assumption of a massless tether removes the effect 
of lateral cable oscillations; and it is under the conditions of low tension (slow rotation 
rates) that were studied that the amplitude of these oscillations can become large. 
Although minimization of the effect of these lateral oscillations could be accomplished 
by making the tether of light material, it is felt that the results presented in this report 
should be considered as qualitative in nature. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 13, 1969. 
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Figure I- Sketch of system. 
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Figure 7 - Angular bound of oscillations of final configuration as a function of Kq where slackness occurred in the initial damped system* 
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(b) k = 2.136 N/ra; E = 2500 m. 
Figure 7.- Continued. 
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Cel k - (XB35 H/m- t l^m 0 m. 
Figure 1- Conclude! 
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